Nonlinear response in overlapping and separated Landau levels of GaAs quantum wells 
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We have studied magnetotransport properties of a high-mobility two-dimensional electron sys- 
tem subject to weak electric fields. At low magnetic field B, the differential resistivity acquires 
a correction Sr oc —\ 2 j 2 /B 2 , where A is the Dingle factor and j is the current density, in agree- 
ment with theoretical predictions. At higher magnetic fields, however, 8r becomes B-independent, 
5r oc — j 2 . While the observed change in behavior can be attributed to a crossover from overlapping 
to separated Landau levels, full understanding of this behavior remains a subject of future theories. 

PACS numbers: 73.43.Qt, 73.63.Hs, 73.21.-b, 73.40.-c 



Among many classes of magnetoresistance 
oscillations 1- - which occur in high Landau levels of 
two-dimensional electron systems (2DES), microwave- 
induced resistance oscillations (MIRO) 1 — are perhaps 
the best known and the most studied phenomenon, both 
theoretically 9 -"— and experimentally^"— In the regime 
of overlapping Landau levels and low microwave power, 
theory predicts that high-order MIRO can be described 
by a radiation-induced correction to the resistivity 
(photoresistivity) of the form 1 ^ 



-Asm2ne UJ , 27re w ^> 1 



(1) 



where p is the resistivity without irradiation, e u = ui/lu c , 
uj = 2nf and lo c are the microwave and cyclotron fre- 
quencies, 



(2) 



A = A X 2 , Ao = 47re u V u — + — 

4t* r 



A = cxp(— 7r/o; c T q ) is the Dingle factor, r q is the 
quantum lifetime, Vu is the dimensionless microwave 
power ; 13 ' 27 r is the transport lifetime, r* is the scat- 
tering time characterizing the correlation properties of 
the disorder potential, 28 and Tj n is the inelastic relax- 
ation time. The first term in the parentheses in Eq. @ 
accounts for the displacement contribution; 9 - " 11 ' 14 ' 29 " — 
owing to the radiation-induced modification of impurity 
scattering, while the second term represents the inelas- 
tic contribution ; 12 ' 14 ' 32 " — originating from the radiation- 
induced change in the electron distribution function. 

Over the past decade, many experiments have ex- 
amined the functional dependences of the MIRO am- 
plitude A on magnetic field i?^ ' 26 ' 36 microwave power 
P^^ILIML™ and temperature T^^. However, di- 
rect quantitative comparison of the measured MIRO am- 
plitude to that predicted by Eq. ^ has not been at- 
tempted to date. The main factor preventing such a 
study is an uncertainty in the microwave power V u ab- 
sorbed by a 2DES. Consequently, it is also not feasible 
to reliably evaluate the scattering parameters entering Aq 
from the measured MIRO amplitude. On the other hand, 



it is indeed very desirable to have a reliable experimen- 
tal probe of such 2DES parameters as and T- m , which 
would allow characterization of the correlation properties 
of the disorder potential and the strength of interactions 
in a 2DES, respectively. 

In this paper we propose and demonstrate an approach 
to experimentally evaluate and r; n in high-mobility 
2DES. More specifically, we employ the nonlinear re- 
sponse of the resistivity to an applied dc field. In con- 
trast to studies investigating the regime of strong elec- 
tric fields ^"-iSii^i^— which is dominated by Hall field- 
induced resistance oscillations (HIRO)j 3 - we focus on the 
regime of weak electric fields. In this regime, to the sec- 
ond order in dc field, the theory^ 4 - predicts, in overlapping 
Landau levels, q dc-induced correction to the differential 
resistivity of the form 



<5r\- 
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(3) 



where Cj = Wj/fou> c , VVj = 2R c e£j is the work done by 
the electric field £j over the cyclotron diameter 2R C , and 



a 



a^X , a — 12ir 



3r 
16tv 



r 



(4) 



Unlike the MIRO amplitude [Eq. @], which contains V u , 
the curvature a [Eq. Q] contains only scattering param- 
eters. 

To examine the applicability of Eqs. ([3]) and (j4]), we 
have measured the differential resistivity in a high- 
mobility 2DES over a wide range of magnetic fields, cov- 
ering the regimes of both overlapping and separated Lan- 
dau levels. At low magnetic fields, we have found that 
the differential resistivity acquires a correction which can 
be well described by Eq. © with a oc A 2 , as prescribed 
by Eq. (TJ|. The obtained value of ao suggests that the 
response is dominated by the inelastic contribution given 
by the second term in Eq. ([J}. At higher magnetic fields, 
we observe a significant deviation from this behavior, 
which we attribute to a crossover from overlapping to 
separated Landau levels. More specifically, at B > 1.3 
kG (1 kG = 0.1 T), we find a = 12tt 2 (B j B f , B « 0.93 
kG. As a result, the correction to differential resistivity 
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FIG. 1. (Color online) Normalized differential resistivity r/p 
versus direct current of density j at different magnetic fields 
from 0.3 kG (top curve) to 1.0 kG (bottom curve) in steps of 
0.1 kG. The maxima of HIRO measured at B — 0.3 kG are 
marked by integers (cf. 1,2). 
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FIG. 2. (Color online) Normalized differential resistivity r/p 
versus ej at different magnetic fields from 0.3 kG to 1.0 kG in 
steps of 0.1 kG (solid lines). Dotted lines are fits to the data 
with r/p = 1 — at 2 j over the range —0.1 < £j < 0.1. 



becomes independent of B and follows Sr/p = —j 2 j j 2 , 
where jo « 4.5 • 1CP 2 A/m. 

Our Hall bar sample (width w — 100 /an) was fab- 
ricated using photolithography from a symmetrically 
doped GaAs/Alo.24Gao.76As 300A-wide quantum well 
grown by molecular beam epitaxy. Ohmic contacts were 
made by evaporating Au/Ge/Ni, followed by rapid ther- 
mal annealing in forming gas. The experiment was per- 
formed in a 3 He cryostat, equipped with a superconduct- 
ing solenoid, at a constant coolant temperature T ~ 1.5 
K. After illumination with visible light, the electron den- 
sity and mobility were n e ss 3.95 x 10 11 cm -2 and 
/i 8.9 x 10 6 cm 2 /Vs, respectively. The longitudinal 
differential resistivity r = dV/dl was recorded using low- 
frequency (a few hertz) lock-in amplification as a function 
of j = I/w at different fixed B ranging between 0.3 and 
2.2 kG. The probing ac current was 0.2 p,A. 

In Fig.Q] we present the differential resistivity r, nor- 
malized to its value at zero current p, as a function of the 
current density j, measured at different magnetic fields 
from B = 0.3 kG (top curve) to B = 1.0 kG (bottom 
curve), in a step of 0.1 kG. The maxima of HIRO, which 
occur at ej = 1,2, are marked by integers (cf. 1,2) next 
to the trace measured at B = 0.3 kG. With increasing 
B, these maxima shift to higher currents and eventually 
move outside the investigated current range. The main 
focus of the present study, however, is the regime of small 
dc fields, €j <C 1, which, according to the theory^ is de- 
scribed by Eqs. ([3]) and ([4]). As seen from Fig.Q] the non- 
linearity in this regime becomes progressively stronger 
with increasing magnetic field. 

Our goal is to analyze the data such as that shown in 
Fig. [1] in terms of Eq. ([3]), extract the curvature a, and 
then discuss it in the context of Eq. (TJ|) . After converting 
the current density j to ej — 2(27r/n e ) 1 / 2 m*j/e 2 B, where 
to* w 0.067toq is the electron effective mass, we replot 



the data shown in Fig.[T] as a function of tj in Fig. [3] 
Presented this way, the differential resistivity shows the 
fundamental HIRO maxima at ej ~ ±1 for all magnetic 
fields, in agreement with previous experimenta l 4 ' 5 ' 42 ' 43 
and theoretica l 13 ' 44 ' 45 studies. Our next step is to fit the 
data with r/p = 1 — ae 2 [cf. Eq. ©] over a range of low 
electric fields, —0.1 < ej < 0.1. Three examples of such 
fits for B = 0.3, 0.5, and 1.0 kG are shown in Fig.rjby 
dotted lines. It is clear that the curvature of the fits, a, 
grows rapidly with increasing B. 

In Fig. [3] we present a parameter a (circles), obtained 
from the fits to the data, such as that shown in Fig. [2 
versus inverse magnetic field 1 / B plotted on a log- linear 
scale. Presented in such a way, the data reveal that the 
parameter a changes by nearly three orders of magnitude 
over the studied B range. The lower B (B < 0.5 kG) 
part of the data can be well described by an exponential 
dependence, a = ao exp(— 27r/w c r q ), in accordance with 
Eq. (TJ|. The slope of our fit to the data (solid line) gener- 
ates T q = 15.2 ps. This value is close to T q obtained from 
the Dingle plots of MIRO and HIRO amplitudes, con- 
firming the validity of our approach. An estimate of ao 
is given by the intercept of the fit with the vertical axis. 
We next analyze the value of ao/127r 2 — 2.25, obtained 
from this intercept, in detail. 

We first recall that the displacement contribution, 
given by 3r + /16r, is sensitive to the correlation prop- 
erties of disorder in the 2DES. For example, for purely 
smooth disorder, the displacement contribution is the 
smallest, t/t* = 12/(r/r q - 1) 0.5 (3t/16t* ps 0.09)4^ 
In the opposite limit of only sharp disorder, t/t* attains 
its maximal possible value, t/t* = 3. We notice that 
even the maximal displacement contribution, 3t/16t* = 
9/16 w 0.56, is small compared to ao/127r 2 = 2.25, 
obtained experimentally. We thus must conclude that, 
regardless of the specifics of the disorder, the inelastic 
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FIG. 3. (Color online) Obtained from the fits (cf. Fig.© of 
parameter a/12n 2 (circles) versus inverse magnetic field 1/B 
plotted on a log-linear scale. The fit (solid line) to the lower 
B part of the data, B < 0.5 kG, with a — Qo exp(— 27r/ai c T q )- 
see Eq. @ -generates r q = 15.2 ps and ao/127r 2 = 2.25. 



contribution dominates the nonlinear response resistiv- 
ity in our high- mobility 2DES. In lower mobility and 
higher density 2DES, the inelastic contribution becomes 
even stronger and the displacement contribution can be 
safely ignored, see, e.g. Refs.|47,48. In our study, the 
ratio Tin/7", which determines the inelastic contribution, 
is bounded by 1.69 < t 1u /t < 2.16, from which we ob- 
tain 0.57 ns < Ti n < 0.73 ns. This result agrees well with 
the theoretical estimate, T m ps 0.56 ns, obtained from 
ft/Tin ^ k 2 B T 2 /E F (E F is the Fermi energy) i 2 . 

We next obtain a more accurate estimate of the dis- 
placement and the inelastic contributions in our 2DES. 
Using t / T s h * 0.2 (t^ is the sharp disorder contribu- 
tion to the quantum scattering rate) obtained from the 
B dependence of the HIRO amplitude^ we estimate 
t/t* ~ 3t/t s „ + 12T q /T ps 0.6 + 0.6 = 1.2, a value 
which reflects approximately equal contributions from 
sharp and smooth components of disorder. 14 Using this 
estimate, we then obtain 3t/16t^ ps 0.23, which leads to 
Tn/T ps 2.0 or Tn ~ 0.64 ns. 

The above analysis shows that the displacement mech- 
anism contributes only a small fraction to the observed 
nonlinearity. Since, theoretically, the relative contribu- 
tions from the displacement and the inelastic mechanisms 
are essentially the same for both the MIRO amplitude, 
Eq. ©, and nonlinear response resistivity, Eq. Q, one 
should indeed expect that the displacement contribution 
to MIRO can also be neglected under similar experimen- 
tal conditions. However, a recent study examining the 
temperature dependence of the MIRO amplitude have 
found no 1 /^-dependence, characteristic of the inelas- 
tic mechanism^ This apparent controversy can be, at 
least partially, resolved by noticing that the density (mo- 
bility) of the 2DES used in Ref.[26| was lower (higher) 
compared to that of the 2DES investigated here. As 
a result, the inelastic contribution, which scales with 



Th/t cx n e / fi, was at least twice as small compared to 
the present study. We also note that under experimental 
conditions of Ref.Hfl the temperature dependence of the 
MIRO amplitude A = X 2 Ao was dominated by an ex- 
ponentially changing A 2 , which significantly complicates 
detecting the temperature dependence of Ti n /r entering 
Aq. To confirm the inelastic contribution in microwave 
photoresistance, it is very desirable to investigate the 
MIRO temperature dependence at lower temperatures 
where A 2 becomes T independent. Such a study, how- 
ever, is complicated by a radiation-induced heating of the 
2DES, which gets progressively stronger at lower temper- 
atures. 

Further examination of Fig. [3] shows that, at higher 
magnetic fields, a grows faster than the exponential de- 
pendence predicted by Eq. (Q} . A departure from the ex- 
ponential behavior is likely a signature of a crossover be- 
tween the regimes of overlapping and separated Landau 
levels. Indeed, using the condition w c T q = 7r/2 ) 50 i 51 we 
find that the Landau levels separate at a magnetic field of 
ss 0.4 kG. Examination of Fig. [3] confirms that the mag- 
netic field, at which the departure from the exponential 
dependence occurs, compares well with this estimate. 

To further examine the regime of separated Landau 
levels, we replot a (circles) in Fig. [4] as a function of B 2 . 
The fit (solid line) to the higher B part of the data, 
B > 1.3 kG, shows that the data in this regime can 
be well described by a/127r 2 = B 2 /B 2 with B ps 0.93 
kG. Since ej cx 1/B, this observation suggests that the 
correction to the differential resistivity becomes B in- 
dependent and is determined only by the applied cur- 
rent j. Indeed, one can write Sr/p — —j 2 /j$, where 
jo = e 2 So V ^/47rv / 67rm* ps 4.5 • 10~ 2 A/m in our 2DES. 

Within a framework of the displacement mechanism, 
such a behavior can be qualitatively understood by not- 
ing that at low electric fields and in separated Landau 
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FIG. 4. (Color online) Obtained from the fits (cf. Fig.© 
parameter q/127t 2 (circles) versus B 2 . The fit (solid line) to 
the higher B part of the data, B > 1.3 kG, with a/127r 2 = 
B 2 /Bl generates B ps 0.93 kG. 
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FIG. 5. (Color online) Normalized differential resistivity r/p 
versus ej measured at B = 0.5 kG at T from 1.5 K to 4.0 K, 
in a step of 0.5 K (solid lines). Dashed lines are fits to the 
data with r/p — 1 — ae 2 over the range —0.1 < ej < 0.1. 



levels, the nonlinear response of the 2DES is governed by 
impurity scattering within a single Landau level, which 
is located near the Fermi surface. In this situation, the 
inter-Landau level spacing is no longer important, and 
the relevant energy scale is given by the Landau level 
width r. As a result, tj — yVj/hu c in Eq. ([3]) should be 
replaced by (3Wj/T, where (3 is a constant of the order of 
unity, which depends on the functional form of the den- 
sity of states. The correction to the differential resistivity 
then takes a form Srj/p oc — W 2 /r 2 . Since Wj oc j and 
is B independent, one obtains Srj/p oc — j 2 /j 2 . There- 
fore, our finding that jo does not depend on B implies 
a B- independent T. While a theoretical expression for 
Sr in separated Landau levels is not currently available, 
we note that Wj G/o) = 2fi.-\/27m e jo/e ~ 0.17 K compares 
well with h/2r q m 0.25 K. 

Finally, we note that the observed nonlinearity weak- 
ens considerably with increasing temperature. In Fig. [5] 
we present normalized differential resistivity r/p versus 
ej measured at B = 0.5 kG at T from 1.5 K to 4.0 K, in 
a step of 0.5 K (solid lines). The fits with r/p = 1 — ae 2 
(dotted lines for 1.5, 3.0, and 4.0 K) demonstrate that Sr 
can still be described by Eq. ([3]) for all temperatures stud- 



ied and that the curvature a decays rapidly with increas- 
ing T. The main source of this decay is the increase of the 
electron-electron scattering, l/rj n oc T 2 . This scattering 
not only suppresses the inelastic contribution, given by 
the second term in Eq. (|4]), but also modifies the quantum 
scattering rate, l/r q i 26 ' ' 52 The latter results in the sup- 
pression of both the displacement and the inelastic con- 
tributions, since both scale with A 2 = exp(— 27r/w c r q ), 
see Eq. . Another source of temperature dependence 
is the enhanced scattering on thermal acoustic phonons, 
which modifies the transport scattering rate, 1/t, and 
gives rise to phonon-induced resistance oscillations j 2 ' 52 i 53 
These oscillations are known to interfere with the nonlin- 
ear response resistivity^^ resulting in a nontrivial, B de- 
pendent corrections to a in Eq. ([3]). Finally, we mention 
a recently reported negative magnetoresistivity effec t 55 1 56 
which occurs in the same range of magnetic fields and is 
strongly temperature dependent. Unfortunately, sepa- 
rating all these contributions does not appear feasible at 
this point. 

In summary, we have studied the nonlinear resistivity 
of a high-mobility 2DES over a range of magnetic fields 
covering the regimes of both overlapping and separated 
Landau levels. At low magnetic fields, we have found that 
the differential resistivity acquires a correction which can 
be well described by Sr oc — exp(— 2ir/ui c Tq)j 2 /lo 2 . Quan- 
titative comparison with existing theory^ indicates that 
the nonlinear response in our 2DES is dominated by the 
inelastic contribution. At higher magnetic fields, we ob- 
serve a significant deviation from the above exponential 
dependence, which we attribute to the crossover from 
overlapping to separated Landau levels. Here, the cor- 
rection to differential resistivity becomes independent of 
B and can be well described by Sr/p = — j 2 /jg, where 
jo w 4.5 • 10~ 2 A/m. It will be interesting to see if future 
theories can explain this finding and clarify the physical 
meaning of B independent jo. 
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